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We present preliminary results for nucleon dipole moments computed with domain wall fermions. Our main 
target is the electric dipole moment of the neutron arising from the 9 term in the gauge part of the QCD lagrangian. 
The calculated magnetic dipole moments of the proton and neutron are in rough accord with experimental values. 



1. Introduction 

One of the most intriguing aspects of QCD is 
that it allows a T and P-odd, gauge invariant in- 
teraction term, the so called 9 term. The pres- 
ence of such a term has the profound effect that 
the Strong interactions violate CP. One of the 
best ways to monitor this is by searching for the 
neutron electric dipole moment, djy. In the Stan- 
dard Model, the CP-odd phase of the CKM mix- 
ing matrix can also trigger a non- vanishing value 
for d,N, but this can not occur at one loop order 
(in the Weak interaction) and is consequently es- 
timated to be < 10 -30 e-cm, many orders of mag- 
nitude below the current experimental boundpQ, 
d,N = |djv| < 6.3 x 10 _26 e-cm. Using this ex- 
perimental bound along with various estimates 
of dpf due to the 9 term then invariably implies 
that the CP-odd parameter in the QCD action, 
9 < 10~ 10 , is exceedingly and indeed unnaturally 
small. Since there is no good reason for this num- 
ber to be so different from unity, its minuteness 
requires extraordinary fine-tuning. This is often 
termed the Strong CP problem. 

To translate the above experimental bound to 
a constraint on the fundamental 9 parameter 
requires evaluation of nucleon matrix elements, 
wherein lattice QCD enters. 

There have been several past attempts to cal- 
culate dj\f , in the continuum and on the lattice. 

1. V. Baluni 2 computed d/v hi the frame- 



work of the MIT bag model obtaining djy 
8.2- 1Q- I6 9e 



cm 



2. Crewther et al. 3 , using an effective chiral 
lagrangian found dN oc 9 M% ln(M^) ~ 5.2- 



10- 16 6»e-cm 



3. Pospelov and Rftz0], using QCD sum rules 
techniques, found dN = 1-2 x 10~ 16 6*e-cm 

4. Aoki and Gocksch|S] were the first to try a 
pioneering lattice calculation of the dN in 
the quenched approximation. 

A non-zero dN is induced by the #-term in the 
QCD gauge action. 
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where Q is the topological charge of the QCD vac- 
uum. Note the factor of i in the 9 term makes the 
action complex and therefore difficult to handle 
in lattice simulations for arbitrary 9. If 9 is small, 
as discussed below, the difficulty can be avoided. 

Using the axial anomaly, one can replace 
the CP violating gauge action above with the 
fermionic action, S' s — —i9m J d 4 xP(x) where 
P{x) — u(x)~f5u(x) + d(x)~f$d(x) + s(x)~f5s(x) and 
to -1 = m^+m^ +m~ 1 . In jS] dN was evaluated 
by extracting the spin-up and spin-down neutron 
masses from the two-point functions obtained by 
adding to the action S' g and a term corresponding 
to a constant background electric field oriented in 



a fixed spatial direction. Since 9 is small, it is suf- 
ficient to consider a single insertion of S' g in the 
2-point functions. 

Since the operator insertion in diagram (a) of 
[5] is manifestly CP-conserving, its contribution 
to a CP-odd physical observable (such as electric 
dipole moment) must vanish. Contributions to 
djv therefore must come from diagram (b) of [H] 
which was ignored in that calculation. 

More detailed reasoning for the fact that cIn re- 
sides only in fig(b) of |o] was given shortly there- 
after by a significant extension of that analy- 
sis appeared recently in 0. As stressed in 0, the 
quenched approximation does not justify ignoring 
diagram (b) of 

2. Our computational strategy 

There are two important elements to our com- 
putational strategy: 

(1) Compute the matrix elements of the 
electromagnetic current between nucleon states, 
(p',s\J^\p, a)g = u(p' 7 s)T^(q 2 )u{p, s), where 



2m 



(2) 



+ (7m 75 Q 2 - 2to75 <7 m ) F A (q 2 ) 
^ » F 3 (q 2 ) 

and use projectors to obtain linear combinations 
of F\ and F 2 , 

Jl - Mq 2 ) + 7^F 2 (q 2 ) (3) 
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and F 3 (q 2 )/2m. By forming ratios of Gj\f(q 2 ) 
and F 3 (q 2 ) with G_g(g 2 ) and taking q 2 — ► 0, 
we find both magnetic and electric dipole mo- 
ments, respectively: e(-Fi(0) + a /1 )/2m and dm = 
eF 3 (0)/2m (note, G E (0) = Fi(0) = 1 and 
for the proton and neutron, respectively, and 
a^ t = F 2 (0) is the anomalous magnetic moment). 
In the above q 2 < 0, and m is the nucleon mass. 

(2) Expand (p' , s\J^\p, s)g to lowest order in 9 
and compute F 3 (q 2 ) in each topological sector v, 
and then average over all sectors with weight Q v . 

A disadvantage of this method is that unlike 
the background-electric-field-method|Sj used in 



[S] , our method does not allow a direct calculation 
of the electric dipole moment, i.e. the value of the 
form-factor at q 2 = since on a finite lattice only 
the form factor F\ in Eq. 2 can be computed at 
q 2 = 9 . Our method requires extrapolation 
of the form factors to q 2 = from non-vanishing 
values of q 2 . 

3. Remarks on the quenched case 

The QCD partition function in the pres- 
ence of explicit CP violation is Z(rj, 77) = 
J VAf, det[D(m) + iOfn^] e^ ^^"^ . Set- 
ting det[D(m) + i9rn r yr > } = 1, we lose 
CP violating physics. However, if 9 is 
small, det [D(m) + i&mys] = det[D(m)] [1 + 
iOrn tr(75D(m) _1 ) ] + 0(9 2 ), and we quench as 
usual by setting det [D(m)] = 1. Considering 
disconnected insertions of P = tr(75.D(m) _1 ), 
it seems clear that one can compute the discon- 
nected diagram (b) of [5] using quenched gauge 
configurations and obtain a non-zero result. 

4. The chiral limit 

The spectral decomposition of ,D -1 (m) leads to 
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for Nf flavors and n + and n_ the number of right - 
and left-handed zero modes of D(m). If we trade 
Q for a disconnected insertion of —mP, <1n will 
vanish in the chiral limit only if (det D(m)) Nf 
vanishes, and thus d/v can not vanish in the 
quenched chiral limit (recall that det D(m) ~ m 
for Q 7^ 0, and contributions to djv vanish for 
Q = 0.). 

5. Numerical Results 

In Table 1 we summarize results for the ratios 
of neutron and proton magnetic to proton elec- 
tric form factors which become the dipole mo- 
ment in question in the limit q 2 — > 0. These re- 
sults were computed on 280 Nf — 2, m sea = 0.02, 
domain wall fermion configurations (separated by 
10-15 trajectories) with m val = 0.04 and 0.08 |TO|. 
The lattice size is 16 3 x 32, L, = 12, and the 



inverse lattice spacing in the m sea — limit is 
a -1 1.7 GeV. We have averaged over time slices 
14-17 and (equivalent) permutations of momenta, 
p= (1, 0, 0), (1,1,0), and (1,1,1). Q was computed 
by integrating the topological charge density after 
APE smearing the gauge fields (20 sweeps with 
ape weight 0.45) We are also investigating 

computing the topological charge from the index 
defined from the domain wall fermion Dirac op- 
erator (strictly valid in the limit L s — > oo). 

Table 1 



Ratios of form factors for m va i = 0.04 (upper) 
and 0.08 ( p proton, N neutron). 



Ratio 


(ap) 2 


value (error) 


G P I /(E + M)/G P 
G P /(E + M)/G P 
G P 4 /(E + M)/G P 


1 

2 
3 


1.363 (57) 
1.392 (61) 
1.540 (91) 


GZ/(E + M)/G P 
G»/(E + M)/G P 
G N M /{E + M)/G P 


1 

2 
3 


-0.836 (38) 
-0.837 (42) 
-0.927(61) 


F»/{E + M)/G P 
F 3 N /(E + M)/G P 
F 3 N /(E + M)/G P 


1 

2 
3 


-0.064 (154) 
-0.037 (157) 
0.150 (172) 


G P M /{E + M)/G P 
G P 4 /(E + M)/G P 
G P M /{E + M)/G P 


1 
2 
3 


1.183 (21) 
1.194 (22) 
1.205 (27) 


GZ/(E + M)/G P 
Gl/{E + M)/G p 
GZ/(E + M)/G p 


1 

2 
3 


-0.739 (15) 
-0.740 (16) 
-0.747 (19) 


F£/{E + M)/G P 
F 3 N /(E + M)/G P 
F 3 N /(E + M)/G P 


1 

2 
3 


-0.048 (53) 
-0.029 (57) 
0.028 (62) 



Similar results hold at m,f = 0.08, but with 
smaller errors. Approximating the q 2 — > value 
of each with the smallest value of p 2 (= (27r/16) 2 ), 
taking rrif = 0.04 as the physical value of the 
light quark mass, and aMjv = aMp = 0.8989(77) 
at this quark mass, we obtain 

a P « 1.45(10) 
< » -1.50(7) 
d N /{ea8) « -0.06(15). 

Considering the crude extrapolations just de- 
scribed, these are roughly consistent with the ex- 



perimental values a P — 1.79 and — —1.91 
(and of course dpj ~ 0).The error estimates are 
statistical uncertainties only. In physical units, 
d N = -7.4(18.0) x lO~ 16 0e-cm, consistent with 
the model calculations mentioned above. We note 
that increases in magnitude for both the pro- 
ton and neutron as m va i decreases. The electric 
dipole moment of the neutron does not show any 
significant dependence on the quark mass, within 
(relatively large) statistical errors. In the near fu- 
ture, we will strive to reduce the statistical error 
on our determination of djy which has already 
yielded an interesting first-principles bound on 
the magnitude of d^. 
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